Vortex States of a Superconducting Film from a Magnetic Dot Array 
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Using Ginzburg-Landau theory, we find novel configurations of vortices in superconducting thin films 
subject to the magnetic field of a magnetic dot array, with dipole moments oriented perpendicular 
to the film. Sufficiently strong magnets cause the formation of vortex-antivortex pairs. In most 
cases, the vortices are confined to dot regions, while the antivortices can form a rich variety of 
lattice states. We propose an experiment in which the perpendicular component of the dot dipole 
moments can be tuned using an in-plane magnetic field. We show that in such an experiment the 
vortex-antivortex pair density shows broad plateaus as a function of the dipole strength. Many of the 
plateaus correspond to vortex configurations which break dot lattice symmetries. In some of these 
states, the vortex cores are strongly distorted. Possible experimental consequences are mentioned. 
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Type II superconductors, with their high critical cur- 
rents and fields, lend themselves to technological applica- 
tions. However, vortices appear in these superconducting 
systems when magnetic fields or currents are made suf- 
ficiently large. Motion of the vortices spoils the perfect 
conductivity important in applications; it is therefore im- 
portant to find ways to pin flux quanta. Systematic stud- 
ies of vortex pinning have been carried out in experiments 
on regular nanoscale arrays where a lattice of defects is 
superimposed on a thin superconducting film. "Antidot" 
arrays, in which pinning centers consist of holes or de- 
pressions in the substrate, were first to be studied [1]. 
Subsequently, magnetic dot arrays have been created, 
usually by the deposition of mesoscopic magnetic dipoles 
on top of the superconducting film [2-7]. In the case of 
the nanoscale magnetic dot arrays, a number of experi- 
mental and theoretical studies have focused on scenarios 
in which each unit cell is penetrated by a finite amount 
of net magnetic flux. In this work, we consider arrays 
of dipoles for which there is no net perpendicular ap- 
plied field. Even without any applied net flux, we show 
that rich vortex phenomena occur as one changes the 
strength of the dot dipoles. We propose an experimen- 
tal scenario to observe these effects, in which an in-plane 
applied magnetic field is used to tilt the dipole moments. 
The thin film geometry prevents the horizontally directed 
field from disturbing the superconducting state, and tilt- 
ing the moments makes possible the adjustment of the 
effective strength of the magnetic dots. 

In our work, we focus on the case in which the dipoles 
prefer an orientation perpendicular to the superconduct- 
ing film. The supercurrents which the magnetic dots in- 
duce move in a clockwise direction. The resulting can- 
cellation of the field from the array magnets is a partial 
realization of the Meissner effect. With the induced cur- 
rents, there is an associated cost in kinetic energy; hence, 
for sufficiently strong dipoles, it is energetically favorable 
to have a vortex in the vicinity of the dot, since this al- 
lows the vortex's counterclockwise currents to partially 
cancel the induced currents. However, due to the zero 



flux condition, vortices cannot form in isolation; vortices 
and antivortices must be generated in pairs. 

Some work has concentrated on isolated magnets, 
where the entire system has cylindrical symmetry [8-10]. 
A few studies in the framework of the London theory 
have focused on a single pair of flux quanta in one unit 
cell, using periodic boundary conditions [11-16]. Depin- 
ning has been studied in the Ginzburg-Landau frame- 
work, but for antidots rather than magnetic dots [17]. 
Because vortices must be put in by hand in the Lon- 
don theory, is it not straightforward to go beyond sim- 
ple situations (generally one vortex-antivortex pair in a 
single unit cell). However vortices appear naturally in 
a Ginzburg-Landau treatment, making tractable more 
complicated situations, such as those involving multiple 
flux quanta pairs in a unit cell, where often one must 
allow for superlattice structures. In this work, to avoid 
assuming the same vortex configuration for each unit cell, 
we study a large (4 x 4) supercell with periodic bound- 
ary conditions. In what follows, p pa ir denotes the num- 
ber of vortex-antivortex pairs per unit cell. We find that 
as the strength of the dot magnets is varied (by tilting 
dipole moments via an in-plane magnetic field) , p pa ir ex- 
hibits sharply defined plateaus. Remarkably, the vortex 
pair density is not always a monotonic function of the 
dipole strength. We will see that there are stable vortex 
states which break orientational and/or mirror symme- 
try. Phase transitions which can be abrupt or continuous 
occur as the dipole strength is varied. In the abrupt case, 
hysteresis phenomena are found, consistent with their be- 
ing first order transitions. Surprisingly, some of these oc- 
cur within a plateau. As will be discussed, second order 
transitions define shifts between plateaus in which p pair 
either increases or decreases via the annihilation or cre- 
ation of a vortex-antivortex pair. A gradual deformation 
of the vortex and antivortex cores is associated with the 
creation or destruction of a flux quantum pair. Since 
London theory does not take into account variations of 
the Cooper pair density, and therefore cannot correctly 
describe vortex cores, Ginzburg-Landau theory is essen- 



1 



tial to describe these novel states. 

Methods and Results- To study our system in the 
Ginzburg-Landau framework, we solve the nonlinear par- 
tial differential equations which one finds on extremizing 
the Ginzburg-Landau internal energy given by 
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In Eq. 1, dimensionless units are used; as a result, all 
linear dimensions are expressed in terms of the super- 
conducting coherence length £. The constant A is the 
condensation energy per unit volume for a uniform bulk 
superconductor and the film thickness (in units of £) is 
given by d. Since d <C £, we view the superconducting 
substrate as a film of negligable thickness. Hence, one 
would not expect an in-plane magnetic field imposed to 
tilt the dipole moments of the dot magnets to affect the 
superconducting state in the film. We assume the meso- 
scopic magnetic dipoles above the substrate to have a 
square cross section. While there is a range of magnet 
thicknesses, we have chosen our dots to be cubic in shape 
with 2.0£ as the length of a side. In our case, the meso- 
scopic magnetic cubes form a square array whose lattice 
constant is 6.25£. We specify the dipole moment of the 
magnets by calculating the positive flux passing through 
each unit cell. Given in units of the fundamental flux 
quantum, this quantity provides a natural measure of 
the dot dipole strength. 

In solving the Ginzburg-Landau equations for our ge- 
ometry, we have replaced continuum variables with their 
discrete versions on a mesh fine enough to ensure conver- 
gence with respect to the discretization (to achieve this, 
we allow at least 5 grid points per coherence length). 
Our scheme of discretization is a gauge theoretic formal- 
ism (discrete versions of the standard continuum gauge 
symmetries are imposed) where currents and vector po- 
tentials Afj and A\- occupy the lattice bonds, while order 
parameter values ipij reside on lattice nodes. We treat the 
mesoscopic magnetic dots as square loops of current with 
a thickness equal to the width of the dot base; the mag- 
netic field and vector potential generated by the magnetic 
dipoles are then easily calculated. Though we handle the 
x and y (in-plane coordinates) discretely, the z direction 
is treated exactly, in the continuum limit. 

Via a conjugate-gradient technique, we solve the 
Ginzburg-Landau equations in an iterative manner. In 
this method, one first linearizes the Ginzburg-Landau 
equations about an initial guess. The solution obtained 
by solving the resulting linear equations is then used as 
an input for the next iteration. We simulate in our calcu- 
lation an experiment in which the effective magnetic dot 
strengths are varied continuously (to realize this in the 
laboratory, one could as noted above apply an in-plane 
magnetic field to tilt the dipole moments of the dot mag- 



nets, thereby varying their effective strength). In "right- 
ward sweeps", we slowly increase the dipole strength. 
Sweeps range from dipoles too weak to generate any 
vortex- ant ivort ex pairs to magnetic dots strong enough 
to destroy the superconductivity altogether. In a right- 
ward strength (toward stronger dipoles), the results of 
each calculation are used as the initial guess of the next 
calculation, in which the Ginzburg-Landau equations are 
solved for slightly stronger magnetic dots. "Leftward" 
sweeps are conducted in a very similar way, with each 
successive calculation using weaker dipoles than the pre- 
vious one. The sweeps in opposite directions complement 
each other by highlighting hysteresis effects, thereby re- 
vealing which transitions are first order. 

To illustrate the configurations which correspond to 
the phases shown in Fig. 2, we display Cooper pair den- 
sities in Fig. 3 and Fig. 4. Vortices and antivortices are 
readily identified as regions of depleted Cooper pairs. Ex- 
ploiting the fact that the order parameter phase is sin- 
gular at the vortex cores makes it possible to distinguish 
vortices from antivortices; the line integral § V(j)-ds yields 
+2-7T (— 2tt) if a (anti) vortex is surrounded by the integra- 
tion contour, and is zero otherwise. In terms of currents 
and magnetic flux, this condition is § j^-ds+&B = ±27r 
for a vortex or antivortex, respectively. This tool for lo- 
cating flux quanta permits the precise depiction of vortex 
configurations, and also conveniently yields ppair- 

Figure 1 is a plot of p pa i r and the Ginzburg-Landau in- 
ternal energy for both sweeps to the right and to the left 
(Egl is given in units of dA^ 3 ). The discrepancy between 
opposite sweeps for p pa ir is readily evident. One can also 
see sudden downward jumps in Eql for both rightward 
and leftward sweeps. These jumps and the hysteresis ef- 
fects in pp a i r are a result of the metastability of some of 
the vortex configurations, signifying that several of the 
transitions are first order. One can construct a "ground 
state energy" by selecting the lowest energy from sweeps 
to the right and to the left, with the preferred state be- 
ing the vortex configuration corresponding to the lowest 
energy. In this manner, we have constructed a phase di- 
agram for the system, shown in Fig 2. Figures 3 and 4 
depict vortex configurations in salient cases. In the phase 
diagram, states are classified according to p pa ir, the num- 
ber of vortex- antivortex pairs. When necessary, as in the 
case of ppair — 2, phases are further classified according 
to the symmetry of the vortex configuration, metastable 
states are not lowest in energy and are therefore not pre- 
ferred by the system. 
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FIG. 1. ppair and Egl plots for sweeps to the left and to 
the right. The energy is given in units of d£ 3 A, where A is 
the condensation energy per unit volume and d is the film 
thickness. 
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FIG. 2. Phase diagram with met ast able states. Stable 
states indicated by numbers in parenthesis; met ast able states 
by numbers in brackets. These numbers indicate how many 
vortex-antivortex pairs per unit cell. Dotted (broken) lines 
indicate energies of metastable states for leftward (rightward) 
sweeps. "Bond" subscripts indicate states in which antivor- 
tices lie along nearest neighbor bonds; "striped" and "diag- 
onal" subscripts indicate antivortices which lie on diagonals 
(next-nearest neighbor bonds). 



antivortex pair per unit cell with vortices in the vicinity 
of the dot centers and antivortices centrally located in the 
interstitial areas. With stronger dipoles, there is a first 
order transition to a configuration for which p pair = 2. 
As can be seen in Fig. 3(a), antivortices (small dark re- 
gions) are connected with the dot magnets by lobes of de- 
pleted Cooper pair density (light areas). These molecule- 
like complexes are aligned along the unit cell diagonals, 
thereby breaking ir/2 orientational symmetry. As the 
dipole strengths are increased further, one finds a state 
which breaks both orientational and mirror symmetries 
[see Fig. 3(b)] . As the sudden symmetry breaking sug- 
gests, the transition is first order. 

Even higher dipole strengths lead to a fractional state, 
for which p pair = 2.5. Figs. 3(b) and 3(c) reveal that 
the transition into this state is a gradual one. Ultimately 
mirror symmetry is restored, as seen in Fig. 3(c). As one 
continues the rightward sweep, one encounters a first or- 
der transition to a p va ir — 3 state [depicted in Fig. 3(d)]. 
Again, in view of the significant dissimilarities between 
the states, this is not surprising. For even larger magnetic 
dot strengths, the vortex pair density exhibits a surpris- 
ing nonmonotonicity by jumping suddenly to zero. De- 
spite the abrupt change in p pa ir, the shift to the p va ir = 
state occurs continuously, through the sequence shown 
in Fig. 4. First, the antivortices rearrange to form the 
Ppair = 3 configuration depicted in Fig. 4(a). Next, the 
antivortices are pulled inward, forming the state shown in 
Fig. 4(b). The triangular structures in Fig. 4(b), each of 
which contains three vortex-antivortex pairs, then begin 
to transform. The vortices move outward from the dots 
to meet the antivortices, which move inward. Ultimately, 
the flux quantum pairs annihilate and the result is the 
Ppair = configuration shown in Fig. 4(c). The final 
transition, one which also occurs in a continuous man- 
ner, is to the pp a i r = 4 configuration shown in Fig. 4(d). 
In each unit cell four vortex-antivortex pairs form along 
nearest neighbor bonds. As the magnetic dipoles are 
made stronger, the pair separations increase, until the 
magnetic dipoles become so strong that the supercon- 
ductivity in the thin film is lost. 

This complicated evolution is driven by the competi- 
tion between intra and interdot potentials for the antivor- 
tices. Evidently, for very strong dipoles, the antivortices 
feel not just the effects of a single magnet, but are affected 
by neighboring magnets as well, and the state deforms in 
such a way that the vortices instead reside on nearest 
neighbor bonds. 



In discussing the phase diagram of Fig. 2, we be- 
gin from the left (weak dipoles) and move to the right 
(toward stronger dipoles). The formation of vortex- 
antivortex pairs is energetically unfavorable for weak 
dipoles and, hence, p pa ir = for the leftmost state. 
The next configuration corresponds to a single vortex- 
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FIG. 3. Cooper pair density plots of stable phases, the im- 
ages in panels (a), (b), (c), and (d) correspond to dot magnet 
strengths of 4.10, 4.62, 4.91, and 5.62 fundamental flux units, 
respectively. Ant ivort ices appear as small dark spots, while 
the large dark spots indicate dot regions; Cooper pair density 
is depressed in regions with lighter shading. 
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FIG. 4. Cooper pair density plots of stable phases. Images 
in panels (a), (b), (c), and (d) correspond to magnetic dot 
strengths equal to 5.85, 5.97, 6.14, and 6.85 fundamental flux 
units, respectively. 

Conclusions and possibilities for experiment- We have 
found, even in the absence of any applied magnetic flux, 
that mesoscopic arrays of magnetic dots can exhibit non- 
trivial vortex phenomena, including configurations which 
break lattice symmetries, states with superlattice struc- 
ture, and a fractional vortex configuration. To realize 
these states in the laboratory, we have proposed an ex- 
periment in which one may vary the effective strength of 
the dot magnets by using an in-plane magnetic field to 
tilt the dipole moment. In this way it should be possi- 
ble to carry out the previously discussed sweeps in mag- 
net strength, making it feasible to access experimentally 
the stable configurations shown in the phase diagram of 
Fig. 2 and even some of the states which we have classi- 
fied as met ast able. 

Finally, we mention some possible experimental tests 
for the novel vortex phenomena discussed in this work. 
Critical currents are a useful probe; different j% and jv 
values would signal broken orientational symmetry. The 
existence of states that break the spatial symmetry of 
the lattice [e.g., Fig. 3(b)] suggests that Ising physics 
might be observed in this system at finite temperatures 
in thermodynamic quantities such as the specific heat. 
It is also interesting to speculate that the non-monotonic 
vortex density found in the vicinity of the four vortex pair 
states (Figs. 4) could lead to a peak in the critical current 
as a function of perpendicular dipole strength. In any 
case, our calculations strongly suggest that this system 
offers a rich phenomenology worthy of further theoretical 
and experimental investigation. 
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